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1.  INTRODUCTION 


This  report  summarizes  work  in  progress  on  the  prob¬ 
lem  of  assimilating  satellite-derived  meteorological  measure¬ 
ments  into  forecast  models  suitable  for  operational  use. 

Tests  have  been  made  with  a  barotropic  and  a  baroclinic  assim¬ 
ilation  model,  both  employing  the  noise-freezing  technique. 

For  the  12-hour  period  between  synoptic  times,  the  assimila¬ 
tion  model  was  run,  taking  in  the  satellite  measurements  at 
the  time  and  place  of  observation.  At  the  end  of  the  12-hour 
period,  the  forecast  of  the  assimilation  model  was  used  to 
alter  the  first  guess  fields  for  the  synoptic  analysis.  In 
this  way  the  information  gained  from  the  satellite  observa¬ 
tions  is  passed  along  to  the  forecast  model,  unless  it  is 
vetoed  by  new  synoptic  data  in  the  analysis  scheme. 

The  new  NEPRF  global  analysis  and  forecast  system  was 
used  to  test  the  barotropic  assimilation  model.  Because  of 
problems  incurred  when  the  new  forecast  system  was  prepared 
for  operational  testing,  and  because  of  limited  computer 
availability,  only  one  forecast  was  made  with  and  without 
the  satellite  data.  The  results  of  this  experiment  were  in¬ 
conclusive,  but  suggest  hope  that  improvement  in  the  forecast 
might  be  demonstrated  by  continued  comparisons. 
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A  baroclinic  assimilation  model  has  also  been  tested, 
but  a  definitive  evaluation  has  not  yet  been  made.  Results 
to  date  indicate  that  the  limited  success  which  was  experienc¬ 
ed  with  the  barotropic  model  may  be  retained  in  the  baroclin¬ 
ic  one. 

Section  2  discusses  the  data  assimilation  problem 
and  the  use  of  the  noise  freezing  technique  in  its  solution. 

The  tests  of  the  barotropic  assimilation  model  with  the  NEPRF 
global  forecast  model  are  detailed  in  section  3,  with  results 
in  section  4.  Some  recommendations  concerning  the  possible 
application  of  the  barotropic  model  are  in  section  5. 

The  baroclinic  assimilation  model  is  described  in 
section  6,  and  section  7  contains  a  discussion  of  some  pre¬ 
liminary  tests. 
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2.  BAROTROPIC  DATA  ASSIMILATION 

The  utilization  of  satellite-derived  temperature 
soundings  in  numerical  prediction  models  is  a  problem  which 
has  defied  solution,  despite  the  efforts  of  many  researchers. 
The  biggest  problem  is  that  the  data  is  distributed  along 
the  time  axis,  and  not  amenable  to  conventional  analysis 
methods.  One  of  the  common  methods  of  assimilating 
satellite  data,  referred  to  as  updating,  involves  replacing 
the  solutions  in  a  forecast  model  with  observed  values  at 
the  time  and  place  of  the  observation.  In  a  primitive 
equation  model,  this  procedure  generates  a  great  deal  of 
gravity  wave  noise  which,  unless  controlled,  can  easily  des¬ 
troy  the  entire  forecast.  If  dissipative  devices  are  used  to 
control  the  noise,  the  quality  of  the  forecast  may  suffer 
more  from  the  dissipation  than  it  is  improved  by  the  data. 

Recently,  the  idea  was  suggested  (Sasaki,  1976b)  that 
the  noise  introduced  by  updating  might  be  controlled  by  al¬ 
tering  the  numerical  model  so  that  the  phase  speeds  of 
gravitational  modes  are  greatly  reduced.  In  a  barotropic 
model,  the  phase  speed  of  gravity  waves  depends  on  the 
square  root  of  the  mean  atmospheric  depth.  By  reducing  the 
modelled  mean  depth,  the  gravity  waves  may  be  slowed  without 
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substantially  affecting  the  meteorologically  significant 


waves.  In  a  baroclinic  model,  the  mean  surface  pressure 
plays  a  similar  role.  This  technique  has  been  referred  to 
as  the  "noise  freezing  method"  (NFM) . 

At  a  particular  synoptic  time,  the  satellite  data  for 
the  previous  12-hour  period  would  be  assimilated  by  running 
the  forecast  model  with  initial  data  from  the  previous  syn¬ 
optic  time,  and  updating  the  model  solutions  with  the  satel¬ 
lite  data.  During  this  assimilation  period,  the  NFM  is 
utilized  to  slow  down  the  noise  waves  and  reduce  their  inter¬ 
action  with  the  large  scale.  At  the  end  of  the  12-hour 
assimilation  period,  the  noise  which  has  accumulated  must 
be  removed,  the  conventional  data  incorporated,  and  the 
forecast  model  is  then  initialized  with  the  resulting  fields 


2.1  MODEL  DESCRIPTION 

A  global  barotropic  primitive  equation  forecast  model 
with  variational  energy  and  mass  conservation  was  the  assim¬ 
ilation  forecast  vehicle  in  the  study.  The  variational  con¬ 
servation  algorithm  is  especially  beneficial  in  the 
assimilation  experiments  to  prevent  instabilities. 

The  model  equations  are: 


a  (hu) 

3  t 


a  cos  0 


D  3 

jy  (huu)  +  j-g-  ( hv v  cos  G) 


This  h  T1  =  0 


hv  (f+H-tan  jj) 

3l 

(i) 
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I  A. 


+  a~ cos ~e  rx  (huv)  +  rs-  (hvv  cos  e)  +  hu  (f+  -^a-  6) 


a  cos  0 


*  i  h  H  -  0 


9  (hu)  ,  9  ,,  . 

—  +  re  (hv  cos  0)  =  0 


where  a  is  the  radius  of  the  earth,  g  is  the  acceleration  of 
gravity,  h  is  depth  of  the  layer,  u  is  the  eastward  wind  com¬ 
ponent,  v  is  the  northward  wind  component,  t  is  time,  X  is 
longitude,  0  is  latitude,  and  f  is  the  coriolis  parameter 
(=  2ft  sin  6,  ft  is  the  angular  velocity  of  the  earth). 


The  model  conserves  the  total  mass  and  the  total 


energy: 


tr  2ir 

/  l 


h  a  cos  9  dX  d0  =  <h> 


TT  2V  ‘ 

/  l  J‘U 


2  +  v2)  +  2-  h2  a2  cos  0  dX  d6  =  <T>(5) 


where  <h>  and  <T>  are  constants  determined  from  the  initial 

fields. 


2.2  VARIATIONAL  CONSERVATION  SCHEME 


A  variational  finite  difference  method  (Sasaki,  1976a) 
is  used  to  solve  the  prognostic  equations  (1)  ~  (3)  under 
the  constraints  (4)  ~  (5).  To  save  computer  time,  an  equal- 
area  longitude-latitude  grid  system  is  used  to  mesh  the 
global  surface.  The  variables  h,  u,  v  are  assigned  at  each 
grid  point. 

The  variationally  designed  scheme  requires  the  follow¬ 
ing  steps:  (i)  One  time  step  forecasting  of  all  the  vari¬ 
ables  by  arbitrarily  selected  finite  difference  form  of  (1) 

-  (3)  (the  Leap-frog-flux  form  is  chosen  in  this  experiment). 
These  forecasts  are  written  h,  u  and  v;  (ii)  Adjusting  h, 
u  and  v  to  satisfy  the  conservation  equations  (4)  ~  (5). 

This  is  done  by  forcing  the  first  variation  of  a  functional 
to  vanish,  i.e. : 


SJ  =  5  I 


5  (u-u) 2  +  a(v-v)2  +  $(h-h) 


+  A, 


I  (^  (u2  +  v2)  +  2.  h2)  -  T° 


+  A 


M 


I  h 


where  6  is  the  first  variational  operator,  J  is  the  function¬ 
al,  I  is  a  sum  over  all  grid  points,  a  and  3  are  constant 
weights,  T°  is  the  sum  of  total  energy  over  all  grid  points 
at  t  =  0,  h°  is  the  sum  of  depth  over  the  entire  grid  at 
t  =  0,  A„  and  A,,  are  the  Lagrange  multipliers  for  total 
energy  and  total  mass  conservation,  respectively.  The 
numerical  method  for  solving  Eq.  (6)  is  described  in 
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Sasaki  (1976a).  Some  numerical  values  are  At  =  5  min., 

A0  =  5°  in  longitude,  A0  =  3.2°  in  latitude  near  the  equator, 
increasing  at  higher  latitudes  in  such  a  way  that  the  area 
of  all  grid  elements  is  constant. 


2.3  NOISE  FREEZING  TECHNIQUE 

Since  the  satellite  data  have  sizable  errors,  con¬ 
siderable  noise  may  result  from  their  use  in  updating  a  fore¬ 
cast  model.  It  seems  highly  desirable  to  prevent  the  noise 
which  results  from  such  errors  from  spreading  so  rapidly. 
However,  we  need  to  ensure  that  the  waves  of  meteorological 
importance  are  not  affected. 

The  basic  idea  of  the  experiment  conducted  to  achieve 
the  above  objectives  will  be  illustrated.  Let  us  consider 
for  simplicity  a  linear  set  of  equations  with  the  assumption 
that  9 /9 y  =  0  except  for  the  coriolis  parameter  f  on  a  local 
cartesian  coordinate  system  (x,y),  which  corresponds  to  the 
complete  form  (1)  ~  (3)  : 


9  u 
9  t 


+  U 


9  u 
9  x 


fv 


(7) 


9  v 
9  t 


+ 


U 


9  v 
9  x 


+  fu 


0 


(8) 


9  h 
9  t 


+  U 


9_h 
9  x 


=  -y  H 


9  u 
9  x 


(9) 


where  U  is  the  mean  zonal  current  velocity,  H  is  the  mean 
depth  of  the  layer  and  y  =  1  (nondimensional  constant) .  The 
corresponding  nondimensional  equations  are 
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_  ,9u'  .  t, ^9  u".  ,  ^ 

R,  (rr-'  +  0  3— >  -  £v  * 


R„  (r^  *v'%]  *  fu'  -  0 


(3H1  +  u"  ^)  = 

3  t'  +  u  9x'; 


where  the  primed  quantities  are  nondimensional  and  the  follow¬ 
ing  quantities  with  asterisks  are  dimensional:  u  =  v*v', 

U  =  v*U',  h  =  h*h',  f  -  f*f',  x  =  x*x',  t  =  t*t';  R  is  the 

o 

Rossby  number  (v*/f*x*)  and  e  =  h*/H.  For  meteorological 
motion,  che  following  relations  are  assumed: 


t*  =  x*/v* 


and  r*  - *  =  1 

f*  x*  v* 


For  global  motion,  R  ranges  from  about  0.1  in  mid-latitudes 

o 

to  about  1.0  in  the  tropics.  However,  e  ~  10  2  since  h*  ~ 
102m,  H  ~  104  m.  Therefore,  the  solutions  are  expanded 
about  the  small  parameter  c : 


h'  =  e°h^  +  e1  h'  +  c2  h'  +  ... 

u  i  2 


u'  =  e°u'  +  e1  u'  +  e2  u'  +  ... 


v'  =  e°v'  +  e1  v'  +  e2  v'  +  ... 


For  e°,  we  obtain  the  same  equations  as  (10)  and  (11)  except 


that  u'  ->  u  ',  v'  -+  v  ',  and  h'  -*■  h  ',  and  (12)  is  simply 
oo  o 

reduced  to  9u'/9x'  =  0.  This  means  that  the  meteorological 
0 

motion  is  primarily  rotational.  The  vorticity  equation  is 
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obtained  by  cross  differentiation  of  (10)  and  (11), 


r  (||; 


+  U' 


9C 


9  x 


-)  -  $'  v'  -  f 


9  u' 
9  x ' 


=  0 


(14) 


where  E,'  =  9v'/9x'  and  0'  =  9f'/9y'.  Elimination  of  9u'/9x 
from  (12)  and  (14)  leads  to  an  equation  equivalent  to  the 
potential  vorticity  equation: 


*  Ro 


+  U' 


9C 


9  h 


t-)  -  +  ef'  + 


U' 


9  h' 
9  x' 


=  0 
(15) 


-2 

Since  e~10  and  y  =  1,  the  vorticity  equation  of  meteoro¬ 
logical  significance  is  obtained  by  omitting  the  terms  mul¬ 
tiplying  e. 


R  (- 


9?. 


9C 


<_  l  ;  ■>  +  U  -~— 

O  9  t  9  x 


^*)  -  e'vo  "  0 


(16) 


Suppose  we  now  arbitrarily  adjust  the  value  of  y.  As  long 
as  y  is  between  1.0  and  0.1,  the  same  result  as  (16)  will 
be  obtained  from  (15).  Since  the  phase  velocity  of  the 
inertial  external  gravity  waves  is  U  +  ZygH  +  f2/k2m  any 
choice  of  a  value  of  y  between  1.0  and  0.1  would  result  in 
slowing  down  the  noise  wave  phase  velocity  without  material¬ 
ly  altering  the  motion  of  meteorological  significance. 

This  statement  must  be  qualified  for  long  waves. 
Reducing  y  in  (9)  decreases  the  effect  of  the  divergence  on 
height  changes.  As  was  shown  by  Cressman  (1959),  the 
correct  specification  of  the  phase  speeds  of  long  waves 
relies  on  the  inclusion  of  the  divergence.  If  y  is  reduced, 
we  may  expect  the  long  waves  to  move  more  slowly. 
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3.  TEST  WITH  NEPRF  GLOBAL  MODEL 

Since  previous  experiments  with  the  barotropic 
noise-freezing  model  produced  better  forecasts  with  the 
satellite  data  than  without  it,  we  decided  to  test  the 
scheme  in  combination  with  the  new  NEPRF  global  model. 
Programs  were  prepared  to  fit  within  the  framework  of  the 
semi-operational  global  model. 

At  each  synoptic  time,  the  barotropic  assimilation 
model  is  initialized  at  the  previous  synoptic  time,  using 
the  500  mb  analyses  for  that  time,  which  are  already  avail¬ 
able.  The  assimilation  model  makes  a  12-hour  forecast, 
absorbing  the  satellite  observations  at  the  proper  time  and 
place.  A  Cressman  weighting  function  was  used  to  spread  the 
influence  of  each  observation  over  a  3-grid  interval  radius. 


where  D  is  the  distance  in  grid  intervals  between  the  grid 
point  and  the  observation. 

Only  the  1000-500  mb  thickness  was  used  from  each 
satellite  sounding.  The  height  of  the  1000  mb  surface, 
interpolated  in  time  and  space  from  the  operational  analyses, 
was  added  to  the  thickness  to  obtain  a  500  mb  height. 


r 
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The  resulting  12-hour  forecast  of  the  500  mb  height 
was  used  to  alter  the  first  guess  fields  for  the  geopotential 
analyses  at  the  next  synoptic  time.  First,  the  500  mb  height 
forecast  was  converted  to  a  1000-500  mb  thickness  by  subtract¬ 
ing  the  1000  mb  height  from  the  first  guess  fields.  For  this 
experiment  these  first  guess  fields  were  the  FNWC  operational 
analyses.  In  an  operational  situation,  they  would  be  the 
12-hour  forecast  of  the  global  NEPRF  model. 

Next,  the  1000-500  mb  thickness  forecasted  by  the 
assimilation  model  was  hydrostatically  converted  to  a  mean 
temperature.  Then  the  1000-500  mb  mean  temperature  of  the 
first  guess  was  computed  and  subtracted  from  the  mean  tem¬ 
perature  forecasted  by  the  assimilation  model.  This  tem¬ 
perature  correction  was  applied  at  every  pressure  level,  and 
corrected  first  guess  geopotentials  were  computed  hydrostat¬ 
ically.  The  temperature  fields  were  not  corrected. 

The  corrected  first  guess  fields  were  then  used  by 
the  analysis  programs  to  generate  the  initial  fields  for  the 
forecast.  Of  course  the  new  synoptic  data  will  overrule  any 
discrepancy  with  the  first  guess. 

So  far  only  one  case  has  been  run  to  evaluate  the 
impact  of  the  satellite  data  assimilation.  Problems  asso¬ 
ciated  with  the  semi-operational  implementation  of  the  NEPRF 
model  have  prevented  further  progress.  The  results  of  this 
test  are  discussed  in  the  next  section. 
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4.  RESULTS  OF  TEST 


The  data  from  1200Z  July  16,  1978  was  selected  at 
random  to  test  the  effect  of  the  assimilation  scheme  on  the 
operational  forecast.  Figure  1  is  the  12-hr.  forecast  made 
by  the  barotropic  assimilation  model,  incorporating  the 
satellite  measurements  as  indicated  in  the  figure.  This  500 
mb  field  was  used  to  alter  the  first  guess  fields  for  the 
1200Z  analysis  as  discussed  above.  Figure  2  is  the  resulting 
500  mb  D-value  analysis  for  1200Z,  and  should  be  compared  to 
Figure  3,  the  analysis  for  the  same  time  without  the  influ¬ 
ence  of  the  satellite  data. 

There  are  several  notable  differences  between  Figures 
2  and  3.  The  most  striking  effect  is  the  low  on  the  east 
coast  of  Argentina  which  is  absent  in  the  analysis  without 
satellite  data.  A  D-value  report  of  -280  m  at  Bahia  Blanca, 
Argentina  was  rejected  by  the  routine  analysis  procedure, 
but  was  accepted  when  the  first  guess  was  altered  by  the 
assimilation  forecast;  a  cyclone  did  develop  in  that  area. 

The  Aleutian  low  was  deepened  by  about  80  m  by  the  influence 
of  the  satellite  data,  while  the  low  off  the  west  coast  of 
Spain  was  filled  by  20  m.  In  the  Southern  Hemisphere  the 
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Figure  1.  12-HR  Assimilation  Forecast  for  12Z  16  Jul78  500  mb  HT  with  satellite  reports 


500  mb  D-value  analysis  for  1200Z  July  16,  78  with  satellite  data  as 


addition  of  the  satellite  data  decreased  the  geopotential 
gradients  south  of  Africa.  The  trough  east  of  Australia 
is  more  pronounced  in  Figure  2. 

The  only  really  striking  difference  between  the  two 
analyses  is  the  low  off  the  east  coast  of  Argentina.  It  is 
to  be  expected  that  this  difference  in  initial  conditions 
should  lead  to  a  large  difference  in  the  NEPRF  model  fore¬ 
casts  made  from  the  two  initial  states.  Unfortunately,  this 
difference  almost  entirely  disappeared  in  the  forecast  model, 
as  we  shall  see. 

Figures  4  and  5  are  the  500  mb  temperature  analyses 
with  and  without  assimilation.  Notice  that  the  Argentine 
low  is  more  prominent  with  assimilation,  but  also  exists 
without  it.  Figures  6  and  7  show  the  500  mb  D-values  after 
variational  balancing.  The  contour  interval  in  these  figures 
is  twice  the  one  in  the  previous  D-value  analyses.  The 
balancing  altered  the  analyses  only  slightly. 

Figures  8  through  31  are  the  forecast  maps  output 
by  the  NEPRF  global  forecast  model.  The  initial  conditions 
are  given  by  Figures  8  through  15.  The  duration  of  the 
forecast  is  given  by  the  first  number  in  the  label.  The 
initial  maps  were  obtained  from  the  balanced  output  of  the 
initialization  programs  by  interpolating  vertically  from  the 
pressure  levels  used  in  the  analyses  to  the  terrain-following 
sigma  levels  used  by  the  forecast  model.  This  interpolation 
involves  a  variational  algorithm  to  maintain  hydrostatic  balance. 
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500  mb  Temperature  Analysis  for  1200Z  July  16,  78  with  Satellite  Data. 


tional  Balancinq  for  1200Z  July  16,  78  with  Assimilation 


Figure  7.  500  mb  D-Value  After  Variational  Balancing  for  1200Z  July  16,  78  Without  Satellite  Data 


With  Assimilation 
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Figure  10.  500  mb  Geopotential  Height  for  S.  Hem.  for  initial  time  from 

NEPRF  Global  Model. 


With  Assimilation 
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Figure  11.  500  mb  Geopotential  Height  for  S.  Hem.  for  initial  time  from 

NEPRF  Global  Model. 
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Figure  16.  24-hr.  Forecast  500  mb  Geopotential  Height  for  N.  Hem, 
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Figure  17.  24-hr.  Forecast  500  mb  Geopotential  Height  for  N.  Hem. 
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Figure  18.  24-hr.  Forecast  500  mb  Geopotential  Height  for  S.  Hem. 
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24-hr.  Forecast  500  mb  Geopotential  Height  for  S.  Hem. 


Without  Assimilation 


Figure  20. 


24-hr.  Forecast  Sea  Level  Pressure  for  N 


With  Assimilation 


■Etf  LEVEL 


24  12Z  1 6\  JUL\ 


Without  Assimilation 


ad  'is/s-tp;  mt-YV  Sf-  L  ■  SEfl  [RlL  PRESSURE— JCPRf /UCLfl.-N 


-V^?  Qi  /,  ' 

m 


h 


*0’**  Ct*T 


Fiqure  30.  48-hr.  Forecast  Sea  Level  Pre 


essure  for  S.  Hem. 


With  Assimilation 


With  Assimilation 


Notice  in  Figure  11  that  after  interpolation  the 


Argentine  low  is  present  at  the  500  mb  level  and  has  prac¬ 
tically  the  same  depth  as  in  the  run  with  satellite  data. 

This  fortuitous  outcome  must  have  resulted  from  the  varia¬ 
tional  algorithm  used  in  the  interpolation.  As  a  result, 
there  is  very  little  difference  between  the  two  forecasts. 
Other  differences  which  appear  in  the  results  of  the  initial¬ 
ization  programs  were  also  reduced  by  the  interpolation. 

This  result  seems  to  suggest  that  the  temperature 
analysis,  especially  at  high  levels,  is  controlling  the  mass 
distribution  which  results  from  the  variational  interpolation. 
Since  the  assimilation  technique  altered  only  the  geopoten¬ 
tials,  its  influence  is  being  overruled  by  the  temperatures 
during  the  interpolation. 

It  can  be  seen  from  Figures  16  through  31  that  the 
two  forecasts  made  with  and  without  the  satellite  data 
influence  are  very  similar.  This  circumstance  obviates  any 
discussion  of  the  accuracy  of  the  forecast,  since  it  is 
obvious  that  no  significant  improvement  was  wrought  by  the 
satellite  data  assimilation. 
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5 .  RECOMMENDATIONS 

The  fact  that  a  significant  improvement  in  the 
fields  resulting  from  the  initialization  model  was  demon¬ 
strated  gives  hope  that  a  positive  impact  on  the  forecast 
might  be  obtained.  Toward  this  end,  we  recommend  that  the 
vertical  interpolation  procedure  in  the  forecast  model  be 
reviewed.  Even  though  the  effect  of  the  interpolation  in 
the  run  without  satellite  data  was  to  improve  the  forecast, 
there  is  no  assurance  that  improvement  will  always  result. 

Further  testing  should  be  done  with  this  scheme  to 
clarify  its  effectiveness. 


6.  BAROCLINIC  PRIMITIVE  EQUATION  MODEL 


A  global  baroclinic  primitive  equation  model  using 
spherical  coordinates  has  been  developed  in  order  to  perform 
the  proposed  research.  This  model  employs  simple  finite  dif¬ 
ferencing  and  uses  variational  algorithms  in  order  to  conserve 
total  energy  and  mass.  Furthermore,  it  has  been  programmed 
using  disk  input  and  output  operations  so  that  only  a  fraction 
of  the  model  occupies  the  central  memory  of  the  computer  at  a 
given  time.  Thus,  it  is  possible  to  run  a  global  model  even 
on  a  computer  with  limited  amounts  of  central  memory.  In  this 
section  the  model  itself  is  described  as  well  as  its  program¬ 
ming  flow. 


Model  Description 

Units:  MGS. 

Vertical  Structure:  KM  specifies  the  number  of  levels  where 

v  and  T  are  carried.  c  =  0  at  250  mb  (P  )  and  o  =  1  at  the 

surface  (P  ) .  o  =  0  at  a  =  0  and  c  =  1.  The  model  consists 
s 

of  KM  layers  each  ^  thick.  The  vertical  structure  is  illus¬ 
trated  in  Fig.  32a. 

Horizontal  Structure:  JM  specifies  the  number  of  grid-points 
around  a  latitude  ring  (e.g.,  JM=48  implies  7.5°  grid  spacing). 
IM  is  set  to  JM/2  +  2  giving  a  constant  latitude-longitude 
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grid  with  two  fictitious  latitude  rings  near  the  poles.  The 
fictitious  rings  are  carried  in  the  model  in  order  to  simplify 
the  coding.  The  Arakawa  (C)  lattice  is  used  as  illustrated  for 
IM=26  and  JM=48  in  Fig.  32b. 

Values  of  <p ,  it  and  a  are  carried  at  the  grid  points 
labeled  T.  The  values  of  v  at  the  north  and  south  poles  are 


seldom  needed  since  cos  0=0  for  0  =  ±90°.  This  will  be  seen 

more  clearly  when  the  primitive  equations  are  written  out  in 

flux  form  for  (1,0, a)  coordinates.  There  are  IM  latitude  rings 

where  T,  <J>,  it,  and  a  are  carried.  For  JM=48,  T.  .  =  T_  . 

1  r  J  I  3  -  ^  ’ 

and  j  =  T25  j  +  24"  T^e  va^ues  °f  <J>  /  7T/  and  a  are  computed 

in  a  similar  fashion  for  the  fictitious  latitude  rings.  The 
values  of  v  at  the  north  and  south  poles,  respectively,  are 


V  —  v 

_  2,i  v2,i±24 


and  v 


25,  j 


-  V24,j  ~  V24,j±24 


Primitive  Equations:  The  primitive  equations  are  written  in 
flux  form  using  (1,0, a)  co-ordinates.  Before  describing  the 
equations  some  preliminary  definitions  must  be  made: 
p  =  iTO  +  pt  where  tt  =  pg  -  pfc, 

T  =  0  (Imo)K'  and 

3  <p 

3F  - 

The  primitive  equations  are: 

Mass  Continuity 

It  +  g-cbs  e  lTX  {7ru)  +  Te  (7r  v  cos  0)1  +  To  *  0  (] 
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Equations  of  Motion 


k  (’nj) 

Tt  "v) 


a 

cos 

0 

TT 

a 

cos 

0 

1 

a 

cos 

0 

[  — 

"v 


ril 

-V 


7T  r3  4>  . 

+  ;  >39  +  °  " 


9  9  • 

( ttu 2 )  +  g-g  (it  u  v  cos  9)]  +  (ima) 

+  CTO  |y]  -  [f  +  H_t|£L_i]TT  V  =  IT  Fx  (18) 

(tt  u  v)  +  g-g-  (it  v2  cos  0)]  +  g-j  (tt  v  a) 

|f]  +  [f  +  ff  U  =  1  Fg  (19) 


Thermodynamic  Equation 


rt  ('T)  +  r^sinr'lx  <™t>  +  re  ^  cos  9)J+  ( 


D  .  <  9  ,  •  . 

rooo’  57  (TOSI 


+  Tr2aa  j- 


1  •r9u,9  ,  ,  3d,  ttO 

- a  t-5-r  +  rs*  (v  cos  0)}  +  5-g]  =  7^-' 

a  cos  09X90  9  a  c 


(20) 


In  the  model  centered-time  and  centered-space  differencing  is 
used  along  with  these  equations  in  order  to  forecast  tt  ,  u,  v, 
and  t. 

The  forecast  for  it  at  time  n+1  must  be  made  first  so 

• 

that  a  diagnostic  equation  for  a  at  time  n  can  be  solved  and 

the  result  used  in  forecasting  u,  v,  and  T  at  time  n+1.  The 

values  of  4>  at  each  time  step  are  found  using  the  hydrostatic 
.9  <J>  . 

equation  («-L  =  -iro)  . 

d  o 

Variational  Conservation  of  Mass  and  Energy:  Using  the  calcu¬ 
lus  of  variations  the  tt,  u,  v,  and  T  fields  are  adjusted  at 
each  time  step  so  that  the  initial  mass  and  energy  input  to 
the  model  are  conserved.  For  energy  the  functional  to  be  mini¬ 
mized  is: 
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J  =  ZEE  5AA  (u-u) 2  +  ZEE  5AA  (v-v)  2  +  ZZZ  BAA  (T-T) 2 


. .  2  •  2 

+  Xg  {ZEE  [c  pirTAA  +  7r3AA^  +  tt1AA^  ]  -  TE0 }  , 

where  TE0  =  ZEECpTr0T0AA  +  EZEtTo^AA^-  +  ZEZtt./aa^- 

(ir0,  T0 ,  u0,  and  v„  are  all  initial  values  and  AA  is  the  area 
increment  associated  with  the  appropriate  grid  point) . 

The  ratio  of  a/B  is  taken  to  be  T/cp  ~  1/3.68.  The 
resulting  corrections  for  u,  v,  and  T  applied  at  each  time 


step  are: 


u  =  u  (1 
v  =  v  (1 


An  approximate  solution  for  A  is  used: 


xe  - 


^  ~  2  •  ~  2 
ZEE(c  ttAAT  +  fr^AA  j  +  tt^A  j  )  -  TEC 


c  2tt2AA  tt^2AAu2  tt12AAv2 


ZEE  (- 


For  mass  the  functional  to  be  minimized  is: 


I  =  SE  AA  ( 7T— tt )  2  +  X  [EE  ttAA  -  EEtt-AA] 

1}  M  ij  lj 

where  tt0  is  the  initial  value  for  tt  .  The  resulting  correction 
to  the  77-field  is: 

XM 

tt  =  tt  -r  ~2~  where 


A  _  2[EETrAA  -  EEtt0AA] 
M  EEAA 


51 


This  correction  is  also  applied  each  time  step.  In 
order  to  apply  the  two  corrections  simultaneously,  is  slight¬ 
ly  adjusted  by  the  ratio  (EETrAA/£l7r0AA)  in  order  to  compensate 
for  the  adjustment  to  tt  and  its  effect  upon  the  energy  calcula¬ 
tions  . 

Filtering:  Two  distinct  and  very  different  types  of  filtering 

are  used  in  the  model.  Robert  time  filtering  is  employed 
to  eliminate  the  "two-solution"  problem  that  leapfrog  time- 
differencing  would  induce.  The  prognostic  quantities  u,  v,  T, 
and  it  are  filtered  each  time  step  using: 

F  =  (1  -  a)  F  +  2.  [p  +  p  ’  ]  and  a  =  0.2. 
n  n  2  n-i  n+l 

In  order  to  ensure  the  computational  stability,  space 
filtering  must  be  performed  in  the  vicinity  of  the  poles  due 
to  the  reduction  in  size  of  the  horizontal  grid  spacing  caused 
by  the  equal  latitude-longitude  grid.  In  this  model  Fourier 
filtering  is  employed  for  each  of  the  prognostic  variables  at 
®^®^y  time  step  for  latitudes  poleward  from  and  including  45°. 
Table  1  indicates  the  maximum  wave  number  passed  for  a  given 
variable  and  latitude  when  IM  =  26,  JM  =  48,  and  At  =  450  sec. 

Although  many  papers  on  the  subject  of  spherical  grids 
and  space  filtering  would  indicate  that  a  larger  time-step  may 
be  used  and  computational  stability  maintained,  it  was  found 
that  this  was  not  the  case  for  a  homogeneous  atmosphere  and 
Phillips'  initial  conditions  for  a  4-wave  global  pattern.  It 
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Table  1.  Maximum  wave  number  passed  during  Fourier  filtering. 


Latitude 


Variable 

86.25 

78.75 

71.25 

63.75 

56.25 

48.75 

2 

6 

10 

14 

18 

21 

wM 

4 

8 

12 

16 

19 

23 

82.50 

75.00 

67.50 

60.00 

52.50 

45 

Latitude 

is  suspected  that  under  less  extreme  conditions  than  those  used 
to  test  this  model  a  larger  time-step  may  be  employed. 

MODEL  PROGRAMMING  FLOW 

In  order  to  make  efficient  use  of  computer  time  and 
space ,  this  model  utilizes  circular  buffers.  Essentially,  only 
six  data  blocks  are  carried  in  core  at  one  time  with  constant 
input  from  and  output  to  disk  being  performed  simultaneously 
and  wherever  possible  covered  by  computation.  Each  data  block 
can  be  thought  of  as  a  one-dimensional  array  dimensioned  by 
NVAR  =  ( 8*KM  +  3)*  JM  and  contains  all  of  the  necessary  variables 
for  one  "latitude  ring"  (actually,  due  to  the  staggered  grid  v' 
is  at  a  different' latitude) .  The  location  within  this  array 
of  a  particular  variable  is  determined  by  the  offsets  described 
in  Table  2. 

The  arrays  AO  and  AN  are  dimensioned  by  4*NVAR  and 
2*NVAR,  respectively.  At  a  qiven  time  array  AO  contains  the 


Table  2 . 


Offsets  of  Variables  in  Arrays  AO  and  AN 


Constant  Field  Levels  Name  of  Offset  Value  of  Offset 


Geopotential  of 
terrain  ht. 

1 

MGZ 

0 

One  Time-level  Fields 

Geopotential  heights, <p 
Vertical  Velocities,  cr 

KM 

KM 

MPHI 

MWV 

JM 

(KM  +  1) 

*  JM 

Two  Time-level  Fields 

Terrain  Pressure, tt  ^ 

1 

MPTNM 

(2*KM 

+ 

1)  *  JM 

Terrain  Pressure, tt 

n 

1 

MPTN 

(2*KM 

+ 

2)  *  JM 

u-wind,  un_^ 

KM 

MUNM 

(2*KM 

+ 

3)  *  JM 

u-wind,  u 

KM 

MUN 

( 3*KM 

+ 

3)  *  JM 

v-wind, 

KM 

MVNM 

(4*KM 

+ 

3)  *JM 

v-wind ,  v 
n 

KM 

MVN 

( 5*KM 

+ 

3)  *  JM 

Temperature,  Tn_^ 

KM 

MTNM 

(6*KM 

+ 

3)  *  JM 

Temperature,  Tfi 

KM 

MTN 

(7*KM 

+ 

3)  *  JM 

three  latitude  rings  needed  to  make  a  forecast  for  the  middle 
ring  while  the  next  latitude  ring  that  will  be  used  is  being 
read  from  disk  into  its  remaining  space.  At  a  given  time  the 
forecast  for  the  middle  latitude  ring  is  being  computed  and 
placed  in  array  AN  while  the  forecast  latitude  ring  for  the 
previous  time-step  is  being  written  to  disk.  The  location  with¬ 
in  arrays  AO  and  AN  of  the  various  latitude  rings  is  also  deter¬ 
mined  by  the  offsets  LMN,  LN,  LPN,  and  LIN  (for  AO)  and  LMNP 
and  LNP  (for  AN) .  In  the  following  combination  of  text  and 
illustrations,  the  programming  flow  of  this  model  will  be  des¬ 
cribed  . 
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At  the  beginning  of  a  particular  time-step  the  disk 

file  IIN  contains  IM  records,  each  record  consisting  of  an  array 

dimensioned  by  NVAR  and  containing  the  previously  described 

information  for  one  latitude  ring.  The  rings  are  carried  in 

order  from  the  north  pole  to  the  south  pole  (starting  and  ending 

with  the  fictitious  rings)  and  will  be  denoted  by  R^ ,  R£»  ...» 

R_w .  Subroutine  BFIRST  reads  the  first  four  rings  into  array 

AO  with  the  offsets  LMN,  LN,  LPN,  and  LIN  assigned  to  R^ ,  R2 > 

R^,  and  R^  ,  respectively.  Using  the  corrections  determined  at 

the  end  of  the  previous  time-step,  the  values  of  tt^,  u^,  v^, 

and  T  in  each  ring  are  adjusted  so  that  total  mass  and  energy 
n 

will  be  conserved.  Then  for  each  ring  the  value  of  is  com¬ 
puted  using  the  hydrostatic  equation.  When  this  has  been 
accomplished  for  each  ring,  fictitious  values  of  v  at  the 


north  pole  are  computed  from  the  values  of  v  in  R^  and  are  in¬ 
serted  in  R2  (see  Fig.  33)  .  Next,  using  the  primitive  equations 
and  R^ ,  R2,  and  R3,  forecasts  for  7rn+1r  un+i'  and  Tn+i  are  made 
for  the  R2  ring  and  are  written  in  the  "n-1"  position  in  array 
AN  with  offset  LNP .  The  forecast  ring  will  be  denoted  by  R2*. 
Thus,  at  this  time  array  R2  contains  the  n-1  and  n  time  values 
of  each  prognostic  variable  while  array  R2*  contains  the  n+1 
and  n  time  values.  Before  the  forecasts  were  made  R2  was  copied 
into  R2*  so  the  forecast  values  (n+1)  simply  replace  the  old 
values  (n-1).  Fourier  filtering  is  now  performed  in  R2*  for  the 
new  forecast  values  (n+1)  of  it,  u,  and  T.  Then,  using  the  n+1 
and  n  values  in  R  *  and  the  n-1  value  in  R, ,  the  n  value  of  R  * 
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Fig. 33.  p.E.  Model  Programming  Flow 

Note:  (u^  denotes  presence  of  7j\,  tv,  <p .  ,  &  ^ . ) 
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is  replaced  by  a  Robert  time-filtered  value.  Computations 
necessary  for  the  determination  of  the  mass  and  energy  constraint 
adjustments  are  now  made  using  the  n.+l  values  in  R2*  and  finally 
the  values  of  Rj*  are  used  to  generate  the  fictitious  ring  R^* 
and  the  ring  R^*  is  written  to  disk  file  IOUT. 

Steps  2  to  IM-3  are  all  similar  in  nature.  Each  step 

begins  with  the  next  latitude  ring  needed  being  read  from  disk 

IIN  into  array  AO,  offset  LIN  and  the  previous  forecast  being 

written  to  disk  IOUT  from  array  AN,  offset  LMNP  (see  Fig.  33)  . 

These  input  and  output  operations  will  be  covered  by  computations. 

Note  that  for  step  2  (and  similarly  for  all  succeeding  steps) 

the  rings  R2,  R^,  and  R^  and  the  ring  R2*  remain  in  the  same 

place  in  arrays  AO  and  AN  that  they  occupied  during  the  previous 

step  while  the  values  of  the  offsets  are  rotated  to  indicate 

their  "new  positions".  The  ring  in  array  AO,  offset  LN  is 

copied  to  array  AN,  offset  LNP  and  forecasts  for  7rn+]_  /  un+i  Vn+1 ' 

and  Tn+^  are  made.  The  appropriate  Fourier  filtering  (if  any) 

is  then  applied  to  the  forecast  values  (n+1)  in  the  R*  ring 

and  the  present  values  (n)  in  the  R*  ring  are  time-filtered 

using  the  n-1,  n,  and  n+1  values  contained  in  the  R  and  R* 

rings.  Using  the  n+1  values  in  the  R*  ring,  the  necessary 

computations  for  the  mass  and  energy  adjustments  to  be  made  on 

the  next  time-step  are  performed.  Finally,  the  mass  and  energy 

adjustments  determined  from  the  previous  time  step  ?re  applied 

to  the  values  of  it  ,  u  ,  v  ,  and  T  in  the  R  ring  which  has  been 

n  n  n  n  3 
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read  into  array  AO,  offset  LIN,  and  the  <J>n  values  for  this  ring 
are  computed  using  the  hydrostatic  equation. 

Subroutine  BLAST  must  be  called  in  order  to  complete 
the  forecast  sequence  for  a  time-step.  First,  -*-s  written 

to  disk  IOUT  and,  using  the  procedure  described  previously, 

rings  RIK_2'  Rim-1'  and  RIM  are  use(*  to  Produce  ring 

(array  AN,  offset  LNP) .  After  this  ring  has  been  Fourier-filtered, 

time-filtered,  and  used  to  compute  mass  and  energy  adjustments, 

it  coo  is  written  to  disk  IOUT.  Finally,  the  values  of  ring  R*IM 

are  determined  from  ring  R*IM_1  and  put  in  array  AO,  offset  LIN. 

From  there  they  are  written  to  disk  IOUT  and  the  sequence  is 

complete. 

Upon  completion  of  the  forecast  sequence  for  a  time- 
step,  the  disk  files  IIN  and  IOUT  are  both  rewound  and  their 
names  switched.  Thus,  the  disk  upon  which  the  forecasts  were 
written  becomes  the  new  input  disk  while  the  old  input  disk 
will  be  used  to  house  the  next  set  of  forecasts.  After  each 
time  step  the  offsets  for  the  two  time-level  fields  are  also 
switched  (i.e.;  MPTNM  and  MPTN,  etc.)  since  during  the  forecast 
process  the  n+1  values  replace  the  n-1  values  while  the  n 
values  remain  in  position.  It  should  be  noted  in  Fig.  33  that 
the  variables  u,  ir ,  and  T  are  carried  on  IM  latitude  rings  (2 
fictitious),  and  v  is  carried  on  IM-1  rings  (including  the 
north  and  south  poles) .  Also,  the  values  of  v  carried  on  a 
given  R  or  R*  ring  are  those  north  of  the  values  for  u,  ir,  and 
T  carried  on  the  same  ring. 
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7.  NOISE  FREEZING  FOR  THE  BAROCLINIC  PRIMITIVE  EQUATION  MODEL 


In  section  2 . 3  it  was  shown  that  the  phase  speed  of  ex¬ 
ternal  gravity  waves  in  a  barotropic  model  could  be  reduced  without 
significantly  altering  the  speed  of  meteorological  waves  when  the 
value  of  y  in  (9)  was  adjusted  from  1.0  to  a  value  between  1.0  and 
0.1.  The  speed  of  the  gravity  waves  is  proportional  to  the  square 
root  of  the  mean  atmospheric  depth,  H.  By  reducing  the  modelled 
mean  depth,  the  gravity  waves  are  slowed  without  substantially  af¬ 
fecting  the  waves  of  meteorological  significance.  Upon  examination 
of  (7) -(9)  we  can  state  that  noise  freezing  is  performed  simply  by 
decreasing  the  contribution  of  horizontal  divergence  to  the  tenden¬ 
cy  of  the  depth  of  the  model  atmosphere,  h. 

The  implementation  of  noise  freezing  in  a  baroclinic 
model  requires  that  the  contribution  of  horizontal  divergence  to 
the  surface  pressure  tendency  be  reduced  without  disrupting  the 
horizontal  gradients  of  n .  This  is  done  most  easily  by  subtract¬ 
ing  a  constant  from  all  values  of  it  carried  within  the  model.  The 
constant  is  chosen  to  be  as  large  as  possible  without  causing  neg¬ 
ative  or  zero  values  of  it.  Since  the  equations  of  motion  (18)  and 
(19)  and  the  thermodynamic  equation  (20)  (repeated  below)  used  in 
the  model  are  written  in  flux  form  by  adding  the  continuity  equation, 
care  must  be  taken  so  that  the  reduction  of  tt  does  not  improperly 
affect  the  values  of  other  variables  in  the  equations.  In  (18) -(20) 


the  variables  indicated  by  arrows  must  not  be  affected  by  the  reduc¬ 
tion  in  7T. 
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When  these  exceptions  have  been  accounted  for,  the  model 
can  be  run  with  noise  freezing  by  simply  subtracting  a  constant 
from  the  initial  values  of  it. 

A  set  of  initial  value  conditions  was  prepared  using 
NMC  analyzed  fields  for  12Z  15  February  1976.  These  data  were 
available  at  all  of  the  standard  levels  at  2  1/2°  intervals  and 
were  interpolated  to  the  grid-points  required  by  our  model. 

Initial  sea-level  pressures  were  estimated  from  the  heights  of  the 
1000  mb  surface.  The  u-  and  v-  components  of  the  wind  were  adjust¬ 
ed  so  that  the  initial  mass  divergence  at  each  grid-point  was  zero 
and  internal  balancing  was  applied  to  adjust  the  geopotential 
(temperature)  fields  so  that  the  initial  mass  divergence  tenden¬ 
cy  was  zero. 
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Twelve-hour  forecasts  were  then  made  with  and  without 
noise-freezing.  The  frozen  forecast  was  made  by  subtracting  675 
mb  from  the  initial  values  of  it  while  accounting  for  this  change 
in  the  determination  of  the  arrowed  variables  in  (18)  —  (20)  .  The 
frozen  and  unfrozen  sea-level  pressure  forecasts  are  shov/n  in 
Figs.  34  and  35,  respectively,  while  Fig.  36  displays  the  verify¬ 
ing  NMC  analysis  for  00Z  16  February  1976  (it  should  be  noted 
that  the  sea-level  pressure  in  Fig.  36  was  also  estimated  from 
1000  mb  heights).  The  largest  difference  observed  in  Figs.  34 
and  35  is  the  intensity  of  the  high  located  at  roughly  55°N,  40°E. 
Both  forecasts  are  quite  reasonable  when  compared  to  the  NMC  anal¬ 
ysis  in  Fig.  36  .  Certainly  the  noise-freezing  does  not  degrade 
the  12-hour  forecast  made  by  the  model. 

Next,  the  model  was  subjected  to  a  large  shock  by  reduc¬ 
ing  the  temperatures  at  all  three  levels  by  20°C  at  33.75°N,  90°E. 
Frozen  and  unfrozen  forecasts  were  run  with  and  without  this  shock. 
Fig.  37  displays  the  difference  field  (shocked  minus  unshocked) 
for  1-hour  unfrozen  forecasts  of  tt  .  The  central  value  of  the  high 
exceeds  48  mb.  The  difference  field  for  3-hour  frozen  forecasts 
of  tt  is  shown  in  Fig.  38  where  the  central  value  of  the  high  is 
just  above  40  mb.  It  can  be  seen  that  noise-freezing  has  reduced 
the  phase  speed  of  the  gravity  waves  by  almost  exactly  a  factor 
of  three.  In  both  figures  a  rather  suspicious  looking  belt  ex¬ 
tends  around  the  globe  centered  approximately  at  45°N.  This  is 
the  direct  result  of  the  Fourier  filtering  performed  in  the  model 
poleward  of  45°N.  A  similar  plot  for  the  3-hour  unfrozen  forecasts 
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12-hour  forecast  sea-level  pressure  field  (mb)  valid  at  002  16 


(mb)  for  00Z  16  Februa 


Fig. 37.  Difference  field  (mb)  between  1-hour  unfrozen  sea-level  pressure  foreca 
(shocked  minus  unshocked)  for  13Z  15  February  1976. 


Fig.  38.  Difference  field  (mb)  between  3-hour  frozen  sea-level  pressure  forecasts 
(shocked  minus  unshocked)  for  15Z  15  February  1976. 


is  displayed  in  Fig.  39  .  It  can  be  seen  that  the  central  value 
of  the  high  has  diminished  to  about  24  mb  and  that  the  radius  of 
influence  of  the  noise  is  roughly  three  times  as  large  as  that 
in  Figs.  4  and  5.  Thus,  it  can  be  seen  in  these  figures  that  the 
noise-freezing  technique  can  be  effectively  applied  to  a  baroclinic 
primitive  equation  model  in  order  to  control  the  noise  induced  by 
assimilation  of  satellite  data  into  the  model  at  asynoptic  times. 

Experiments  are  in  progress  to  determine  whether  the 
baroclinic  assimilation  model  can  produce  initial  fields  which 
result  in  an  improved  forecast.  We  are  using  the  data  which  was 
collected  during  Data  Systems  Tests  5  and  6  during  the  summer  of 
1975  and  the  winter  of  1976.  These  data  sets  provide  unprecedented 
coverage  of  satellite  and  other  asynoptic  observations.  Figure  40 
depicts  the  global  coverage  by  satellite  soundings  for  a  12-hour 
period.  Figure  41  shows  positions  of  satellite  wind  measurements 
made  by  NMC .  For  these  test  periods  a  group  at  the  University  of 
Wisconsin  made  wind  measurements  from  satellite  pictures  using  a 
special  technique.  Their  coverage  is  in  figure  42  . 

The  wind  data  presents  a  special  problem  to  the  assimi¬ 
lation  scheme  because  all  the  measurements  for  the  12-hour  period 
were  made  at  only  two  times. 

Forecasts  with  and  without  the  satellite  soundings 
shown  in  figure  40  were  made  for  the  12-hour  period.  Figure  35 
is  the  sea-level  pressure  forecast  without  assimilation.  Figure 
43  shows  the  assimilation  forecast  without  the  use  of  noise  freez¬ 
ing,  and  figure  44  is  the  equivalent  forecast  with  the  satellite 


Difference  field  (mb)  between  3-hour  unfrozen  sea-level  pressure  forecasts 
(shocked  minus  unshocked)  for  15Z  15  February  1976. 
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Locations  of  SIRS  and  NIMBUS  6  Satellite  Soundings  for  the  12-hour  period  1200Z 


Positions  of  NMC  Satellite  Wind  Measurements  for  the  12-hour  period  1200Z  to 
2400Z,  15  Feb.  ,  1976.' 


Fig.  42.  Position  of  Satellite  Wind  Measurements  made  by  University  of  Wisconsin  for  the 
12-hour  period  1200Z  to  2400Z,  15  Feb.,  1976. 


data  and  noise  freezing.  There  are  significant  differences  in 
these  forecasts,  particularly  in  the  Southern  Hemisphere  and  near 
Alaska.  The  difference  fields  showing  the  influence  of  the  satel¬ 
lite  soundings  are  shown  in  figures  45  (without  noise  freezing) 
and  46  (with  noise  freezing) . 

It  seems  that  the  most  important  role  of  the  noise 
freezing  is  to  slow  the  geostrophic  adjustment  mechanism,  so  that 
the  influence  of  the  asynoptic  data  is  retained  during  the  12-hour 
period.  No  attempt  has  yet  been  made  to  evaluate  the  accuracy  of 
these  forecasts. 
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